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ABSTRACT 

The Gaia-ESO Survey has recently unveiled the complex kinematic signature of the Gamma Velorum cluster: this cluster 
is composed of two kinematically distinct populations (hereafter, population A and B), showing two different velocity 
dispersions and a relative ~ 2 km s“^ radial velocity (RV) shift. In this paper, we propose that the two populations of 
the Gamma Velorum cluster originate from two different sub-clusters, born from the same parent molecular cloud. We 
investigate this possibility by means of direct-summation N-body simulations. Our scenario is able to reproduce not 
only the RV shift and the different velocity dispersions, but also the different centroid (~ 0.5 pc), the different spatial 
concentration and the different line-of-sight distance (~ 5 pc) of the two populations. The observed 1 — 2 Myr age 
difference between the two populations is also naturally explained by our scenario, in which the two sub-clusters formed 
in two slightly different star formation episodes. Our simulations suggest that population B is strongly supervirial, while 
population A is close to virial equilibrium. We discuss the implications of our models for the formation of young star 
clusters and OB associations in the Milky Way. 

Key words. Methods: numerical - Stars: kinematics and dynamics - Galaxy: open clusters and associations: individual: 
the Gamma Velorum cluster - Stars: formation 


1. Introduction 

The ongoing Gaia-ESO Survey (GES) at the VLT is a gold¬ 
mine of knowledge concerning the formation and evolution 
of young star clusters (YSCs) and OB associations. The 
aim of the GES, which began on December 31 2011 and 
will be completed in ^ 5 years, is to obtain high quality, 
unifor mly calibrated spectra of > 10® stars in the Milk y 
Way (|Gilmore et al.l 120121 : iRandich. Gilmore et al.l 120131) . 
The GES targets include a significant number of YSCs. 
The radial yelocity and the chemical composition is being 
deriyed for a large number of YSC members with high ac¬ 
curacy. This info rmation, when combin ed with data from 
the Gaia mission (jPerryman et al.l|200lh . will enable recon¬ 
struction of the full three-dimensional spatial distribution 
and kinematics of seyeral YSCs. 


The first GES YSC target was the Gamma Velorum 
cluster, a marginally bound cluster of low-mass pre-main 
sequence (PMS) stars, surrounding the massiye Wolf-Rayet 
(WR) - O binary sy s tem named Velorum (HD 68273 , 
WRll; ISmithI 119681 : ISchaerer. Schmutz fc Greii^ I1997D . 
The Wolf-Rayet and O star components of 7 ^ Velorum 
haye current masses ^ 9 ± 2 and 30 ± 2 M 0 , respectiyely 
(iDe Marco fc Schmutd . [19991) . which impli es that the initia l 
masses were 35 and 31 Mq, respectively (lEldri d^ .[ 200 ^ . 
7 ^ Velorum is the most massive member of the common 
proper motion Vela OB2 association, comp osed of 93 early- 
type stars, spread over 100 square degrees (|de Zeeuw et al.l 
19991). The Gamma Velorum cluster was firstly identified by 
Pozzo et (l2000r) . on the basis of the strong X-ray emis- 
sion of the PMS s tars, and then further investigated by 
I Jeffries et al.l (|2009l) . The distance of the Gamma Velorum 
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cluster is about 350 pc (iPozzo et al. 120001: 

van LeeuwenI 

120071: iMillonr et al.l l2007l North et al.l 12007 

). The PMS 


stars in the Gamma Velorum cluster have an estimated age 
~ 5 — 10 Myr (or even more), which is only marginally con- 
sistent with th e age of the 7 ^ Veloru m binary 5.5 Myr, 
lEldrideel 1200911 ■ iJeffries et ahl (l2014ll suggest, on the basis 
of the Lithium (Li) depletion observed among its M-dwarf 
members, that the Gamma Velorum cluster is older than 
10 Myr. 

The GES data (|Jeffries et al.1 l2014ll show that the 
Gamma Velorum cluster has a complex kinematic struc¬ 
ture: it is composed of two kinematic ally distinct popula - 
tions (named population A and B bv IJeffries et al.l[2014ll . 
Population A (hereafter pop. A) has a 2.15 ± 0.48 km s“^ 
lower radial velocity (RV) than population B (hereafter 
pop. B). The velocity dispersion of the two populations is 
also different: pop. A (pop. B) has a best-fitting velocity 
dispersion cta = 0.34 ± 0.16 km s“^ (ctb = 1-60 ± 0.37 km 
s-i). 

The colour-magnitude diagram suggests that either 
pop. A is younger than pop. B by 0.4 ± 0.6 Myr (if the two 
populations are at the same distance), or pop. A and pop. B 
are coeval but pop. A is closer to us by 4±5 pc. The spectral 
analysis indicates a possible age difference, but in the oppo¬ 
site sense: pop. A appears to be marginally more Li depleted 
than pop. B, which implies that pop. A might be 1 — 2 Myr 
older. The age difference derived from Li depletion, com¬ 
bined with the analysis of the colour-magnitude diagram, 
hints that pop. A might be closer to us than pop. B by a 
few parsecs. 

Finally, the spatial distribution also suggests differences 
between pop. A and B: the former is slight ly more con¬ 
centr ated than the latter (see Figure 11 of IJeffries et al.1 
I 2 OI 4 I) . and the centroids of the two populations are offset 
by 4.4 ±3.0 arcmin (0.4 ±0.3 pc for an assumed distance of 
350 pc). The centroid of pop. A is almost coincident with 
the position of 7 ^ Velorum, indicating a connection between 
po p. A and the WR-0 binary system. 

IJeffries et al.l (|2014ll propose three possible scenarios for 
the formation of the Gamma Velorum cluster. In the first 
scenario, pop. A and pop. B are the remnants of a larger 
cluster. After the evaporation of the parent molecular gas, 
only pop. A remained bound and in virial equilibrium, while 
pop. B became unbound. This scenario fails to explain the 
difference in RV between the two populations and the differ¬ 
ent Li depletion. In the second scenario, 7 ^ Velorum formed 
in a supervirial OB association. Pop. A represents the com¬ 
ponent of the initial OB associati on which was gravi tation- 
ally captured by 7 ^ Velorum ('e.g. lParker et al.ll2014ll . while 
pop. B remained unbound. This second scenario explains 
the velocity differences, but cannot easily explain the age 
differences betwee n the two populatio ns. In the third sce¬ 
nario, favoured bv iJeffries et al.1 (l20I4^ . pop. A formed in a 
denser region, while pop. B formed in a less dense region: 
pop. A is still marginally bound, while pop. B expands in 
the Vela OB2 association. This scenario could explain an 
age difference between the two populations, as well as a dif¬ 
ferent spatial position and a different RV, because pop. A 
and pop. B formed in two slightly different star formation 
episodes. 

In this paper, we focus on the third aforementioned 
scenario, investigating its implications in detail, by means 
of an N-body model of the Gamma Velorum cluster. 
The main idea is the following: simulations of turbulent 


molecular clouds suggest that star clusters form fr(3m the 
merger of sub-clusters (iBonnell et al.1 20031: Ba^ 20091 
Bonnell et al.l 120111 : iGirichidis et al.1 2011 : Kruiissen et al.l 


2012fl . Furthermore, observations of star clusters embedded 


in molecular clouds indicate that the process of star forma¬ 
tion i s highly sub-structu r ed (e.g . iGartwrigh t fc Whitworth! 
2004 iGutermuth et al.l 2009^ ISanchez fc Alfarol 120091 : 


Andre et alll2010l: iDe Marc hi et al.1120131: iGouliermis et alJ 


2014 : Kuhn et al.l 1201^ Kirk. Offner fc Redmond 2014 ). 

The merger of sub-clusters can also explain some intrigu¬ 
ing dyna mical properties of YSGs (e.g. a fast mass seg¬ 
regat ion, iParker et al.ll20l3: iFuiii et ^20121 : iParker et al.l 
120141) . Thus, we propose that the two populations of the 
Gamma Velorum cluster originate from two different sub¬ 
clusters, born in the same parent molecular cloud. The ve¬ 
locity shift between two sub-clusters is expected to be of the 
same order of magnitude as the turbulent velocity in the 
parent molecular cloud. The observed turbulent motions of 
Galactic molecular clouds are of t he order of ^ 1 — 2 km 
s“^ (e.g. iMcKee fc Ostrikeill2007l and references therein), 
significantly larger than the sound speed of cold gas (~ 0.2 
km s~^). Furtherm ore, stars in the Orion Nebula Gluster 
(|Tobin et al.ll2009l) and molecular gas i n star-forming re¬ 
gions such as Perseus and Ophiuchus (|Andre et all 120071 : 
iKirk et aI1l2007l:lRosolowskv et al.ll2008D exhibit large-scale 
linear velocity gradients ( ;$, 1 km s“^ pc“^). Velocity gra¬ 
dients (~ 0.2 — 2 km s“^ pc“^) have also been found in 
simulations of mo lecular clouds and YSC formation (e.g. 
lOffner et al.1 l2009l) . Such velocity offsets would naturally 
explain the observed RV difference of 2 km s“^ between 
pop. A and pop. B. We investigate this possibility by means 
of direct-summation N-body simulations. 

This paper is organized as follows. In Section 2, we dis¬ 
cuss the simulation method. In Section 3, we present our 
results. In Section 4, we discuss the main advantages and 
drawbacks of our simulation method, and we propose future 
simulations and observational tests. 


2. Methods and simulations 

In this paper, we simulate the collision between two 
different sub-clusters, by means of direct-summation 
N-body simulations. The simulations were performed 
with the STARLAB publ ic software environment 
(iPortegies Zwart et al.L l200lD . in the version modifie d 
by Mapelli et al. (2013; see also iMaoelli fc Bressanl[2013[) . 
STARLAB integrates the dynamical evolution of a star 
cluster, by means of a predictor-corrector Hermite scheme 
(implemented in the kira routine), and follows the 
stellar and binary evolution of the star cluster members 
(implemented in SEBA). 

Initially, the two sub-cluster s are described as two 
Plummer spheres (lPlummerlll91ll) . The single stars and the 
primary members of binaries were randomly generated fol- 
lowing a Kroupa initial mass function (IMF, iKroupalbOOll) 
with minimum and maximum mass 0.1 and 150 Mq, respec¬ 
tively. The secondary members of the binaries were ran¬ 
domly generated from a uniform distribution between 0.1 
Mq and TOi (where mi is the mass of the primary mem¬ 
ber). We enable stellar and binary evoluti on at solar metal- 
licity (consistent w ith Spina et ^ 12014) . as described in 
IMaoelli fc Bressaiil (l2013l) . In the following, we call pop. A 
and pop. B the population of the first and the second sim¬ 
ulated YSGs, respectively. 
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Table 1. Initial conditions of the body simulations. 


Run 

Na 

Ma (Mq) 

rA (pc) 

Qa 

Nb 

Mb (Mq) 

Tb (pc) 

Qb 

/bin 

At! (km s 

D (pc) 

mi, m 2 (Mq) 

Run 1 

780 

506 

2 

1.0 

585 

271 

1 

4.5 

0.46 

2.4 

5 

32, 14 

Run 2 

780 

510 

2 

0.5 

585 

271 

1 

4.5 

0.46 

2.4 

5 

33, 16 

Run 3 

780 

514 

2 

2.0 

585 

271 

1 

4.5 

0.46 

2.4 

5 

39, 6 

Run 4 

780 

513 

2 

3.0 

585 

271 

1 

4.5 

0.46 

2.4 

5 

15, 15 

Run 5 

780 

513 

2 

1.0 

585 

276 

1 

0.5 

0.46 

2.4 

5 

32, 14 

Run 6 

780 

513 

2 

1.0 

585 

277 

1 

2.0 

0.46 

2.4 

5 

32, 14 

Run 7 

780 

513 

2 

1.0 

585 

273 

1 

12.5 

0.46 

2.4 

5 

32, 14 

Run 8 

780 

506 

2 

1.0 

585 

271 

1 

4.5 

0.46 

1.5 

5 

32, 14 

Run 9 

780 

506 

2 

1.0 

585 

271 

1 

4.5 

0.46 

2.0 

5 

32, 14 

Run 10 

780 

506 

2 

1.0 

585 

271 

1 

4.5 

0.46 

3.0 

5 

32, 14 

Run 11 

780 

506 

2 

1.0 

585 

271 

1 

4.5 

0.46 

3.5 

5 

32, 14 

Run 12 

156 

99 

2 

1.0 

117 

56 

1 

4.5 

0.46 

2.4 

5 

5, 4 

Run 13 

156 

97 

2 

4.0 

117 

48 

1 

8.0 

0.46 

2.2 

5 

3, 3 

Run 14 

780 

506 

2 

1.0 

585 

271 

1 

4.5 

0.46 

2.4 

10 

32, 14 

Run 15 

780 

512 

2 

1.0 

585 

295 

1 

4.5 

0.0 

2.4 

5 

— 

Run 16 

800 

512 

2 

1.0 

560 

271 

1 

4.5 

0.75 

2.4 

5 

57, 32 

Run 17 

780 

491 

1 

1.0 

585 

271 

1 

4.5 

0.46 

2.4 

5 

8, 1 

Run 18 

780 

506 

2 

1.0 

585 

272 

2 

4.5 

0.46 

2.4 

5 

32, 14 

Run 19 

1040 

677 

2 

1.0 

585 

271 

1 

4.5 

0.46 

2.4 

5 

56, 35 

Run 20 

780 

506 

2 

1.0 

780 

359 

1 

4.5 

0.46 

2.4 

5 

32, 14 

Run 21 

780 

506 

2 

1.0 

585 

271 

1 

4.5 

0.46 

-2.4 

5 

32, 14 


Run (column 1): identifying name of the run; Na (column 2): number of particles in the first cluster (identified with pop. A); Ma 
(column 3): initial mass of the first cluster; ta (column 4): scale radius of the first cluster; Qa (column 5): virial ratio of pop. A; 
Nb (column 6): number of particles in the second cluster (identified with pop. B); Mb (column 7): initial mass of the second 
cluster; tb (column 8): scale radius of the second cluster; Qb (column 9): virial ratio of pop. B; An (column 10): initial velocity 
shift between the centres of mass of the two clusters; D (column 11): initial distance between the centres of mass of the two 
clusters; mi, m 2 (column 12): initial mass of the members of the most massive binary system in each run. The 21 runs listed in 
this Table are a selected sample among all the runs we performed. They represent the runs that best match the observations 
(runs 1, 2, 13 and 21), and/or that are important to understand the influence of various parameters on the results. 



-14 -7 0 7 14 -14 -7 0 7 14 

X (pc) X (pc) 


Fig. 1. Position of the simulated stars in the xy plane at 
t = 0, 2.7, 4.8 and 6.9 Myr in run 1. Red crosses: cluster 1 
(corresponding to population A), black open circles: cluster 
2 (corresponding to population B). The centre of the frame 
is the centre of mass of cluster 1 (pop. A). 


We ran several simulations with different initial condi¬ 
tions. The most relevant runs and the corresponding pa¬ 


rameters are summarized in Table [T] The most important 
initial parameters appear to be (i) the number of particles 
in each star cluster [Na and iVe, respectively) and the cor¬ 
responding initial masses (Ma and Mb, respectively); (ii) 
the Plummer scale radius of each star cluster [ta and rp, 
respectively); (iii) the relative velocity (Au) and the initial 
distance {D) between the centres-of-mass of the two clus¬ 
ters; (iv) the binary fraction /bin; (v) the deviation of the 
star cluster from the virial equilibrium. The last condition 
is particularly important, because several o bserved YSCs 
show indication of non-virial velocities (e.g. ICottaar et al.l 
I 2 OI 2 D . To quantify the deviation from virial equilibrium, we 
define the Q parameter asQ = K/\W\ (where K and W are 
the kinetic and the potential energy of the system, respec¬ 
tively). For a system in virial equilibrium Q = 0.5. Systems 
with Q > 0.5 (i.e. supervirial systems) might be the out¬ 
come of the evaporation of the parent gas. In fact, the gas 
component is not included in our models, but simulating 
a stellar population with initial velocities above the virial 
value mimics the fact that the potential well changed very 
fast in the recent past, i.e. that the parent gas evaporated 
rapidly. 

We test two extreme values for the initial mass: a total 
mass (pop. A-bpop. B) of ~ 150 M 0 (runs 12 and 13) and 
of ~ 800 Mq (the other runs). The mass of the clusters 
simulated in runs 12 and 13 is similar to the current to¬ 
tal mass of the Gamma Velorum cluster (where for current 
total mass we mean the mass within the 0.9 deg^ area sur¬ 
veyed by GES), while a total mass of ~ 800 M 0 is likely 
closer to the initial mass of the system. The last column of 
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Table [T] shows the most massive primordial binary system 
in each run. A binary as massive as 7 ^ Velorum is likely 
to form only in the most massive systems amongst those 
simulated here (~ 800 Mq). The large scatter in the val¬ 
ues of Ma and Mb listed in Table [T] (given a fixed value of 
Na and Nb, respectively) is due to stochastic fluctuations, 
because we randomly generate the masses of star particles 
according to a Kroupa IMF. 

In run 21 the two sub-clusters do not collide. They are 
initially located at ZI = 5 pc (Table [J), and their separa¬ 
tion increases (negative Av). This run has been performed 
to check whether the properties of the Gamma Velorum 
cluster are consistent with a mere superposition between 
two sub-clusters (without physical interaction). 


3. Results 

We simulate the collision between two different YSCs, 
during and after the first approach. Figure [T] shows the 
time evolution of run 1 in the xy plane. In all runs, the 
initial velocity shift is entirely along the y axis. In run 1, 
the two clusters are initially at a distance of 5 pc along 
the y axis (which is less than the typic al radius of a gian t 
molecular cloud in the Milky Way, e.g. iDame et al.ll2001^ . 
They collide at time t ~ 2.5 Myr. The YSCs simulated in 
the other runs have similar evolution, but the collision time 
is different, since it depends on the relative initial distance, 
velocity and masses. 

Detailed information about all runs (including those 
shown in Figures [T] and [2]) can be found in Table [2j In 
calculating the RVs in Tabled] (as well as in Figured]), 
we have assumed that the line-of-sight is along the y axis. 
Thus, the RV is the velocity along the y axis, i.e. the axis 
along which we generated the velocity shift Av between the 
centre-of-mass of the two clustertQ. Therefore, the plane of 
the sky is the xz plane of the simulations. 

In Table dl we consider as reference times t ^ 2.5, 5 
and 7 Myr since the beginning of the simulations. We stress 
that we select these three reference snapshots, since the ra¬ 
dial distribution of the simulated stars, and especially the 
shift along the line-of-sight between pop. A and pop. B, 
are not consistent with the observations if t << 2.5 Myr or 
t » 7.5 Myr (see Section 15^ . On the other hand, the RV 
distribution does not change significantly at different times 
(see Tabled]), provided that t > 2.5 Myr (i.e. provided that 
the two sub-clusters already collided). Thus, the main con¬ 
straints on the time elapsed since the collision come from 
the radial distribution, while the RV distribution is less af¬ 
fected. Finally, we recall that the time indicated in Tabled] 
does not necessarily correspond to the age of the stars in 
pop. A and pop. B: the tabulated time is the time elapsed 
since the beginning of the simulations, but the two sub¬ 
clusters might have formed before t = 0. Evaluating the 
age of the Gamma Velorum cluster, which is rather un¬ 
certain (see the Introduction), is beyond the aims of this 
paper. 


^ Without this requirement, the probability of observing the 
two populations nearly superimposed would be much smaller, 
even if this is would be compensated by a less stringent require¬ 
ment on the line-of-sight direction. See Section 13.41 for details. 


3.1. Radial velocity and proper motions 

Figured] shows the distribution of the RVs for nine selected 
runs at t = 4.8 Myr (for other runs, see Table dj)- The 
simulated RVs shown in Figured] include observational un¬ 
certainties. These were random ly drawn from the distri- 
bution reported in Figure 2 of l.leffries et al.l (1201411 . The 
zero-point of the simulated RVs is the RV velocity of the 
centre-of-mass of the total system (pop. A+pop. B stars). 
In Figured] we also s how the observed R V distribution (the 
same as in Figure 6 of iJeffries et al.l[2014ll . The observed RV 
distribution has been shifted to match the simulated ones. 

In our simulations, we exactly know which stars belong 
to the first and to the second cluster, respectively, and we 
can derive their intrinsic velocity dispersion^ Columns 7 
and 8 of Table [2] show the intrinsic velocity dispersions (oa 
and ctb) of the simulated pop. A and pop. B, along the line- 
of-sight. Column 9 of Table d] (ARV) shows the intrinsic 
velocity shift of the two simulated populations along the 
line-of-sight, i.e. the difference between the mean velocity 
of pop. A and that of pop. B along the line-of-sight. 

On the other hand, we want to fit the to tal RV distri¬ 
bution of each simulation in the same way as iJeffries et al.l 
(|2014fl . for a better comparison with the data of the Gamma 
Velorum cluster. Thus, we apply a fit with two Gaussian 
components to the total pop. A+pop. B distribution, for all 
runs. The fitting procedure is carried out on binned data, 
and adopting a procedure equivalent to the one described in 
section 4 of Jeffries et al. (2014, i.e. statistically correcting 
for the binary fraction and for the observational uncertain¬ 
ties a posteriori). From this fitting procedure, we obtain the 
velocity dispersions of the two simulated populations along 
the line-of-sight {a a and ctb), the velocity shift of the two 
simulated populations along the line-of-sight (ARV), and 
the fraction of estimated pop. A members with respect to 
the total number of simulated stars (/a)- The results are 
shown in columns 3, 4, 5 and 6 of Table [5] The intrinsic 
velocity dispersions generally agree with the ones derived 
from this fitting procedure, but they might also disagree by 
a factor of ^ 2 (especially for pop. B), when the distribu¬ 
tion of simulated RVs deviates too much from a Gaussian 
distribution, and when it is more difficult to disentangle 
pop. A from pop. B. 

We adopt two criteria to check whether a simulation is 
consistent with the observations of the Gamma Velorum 
cluster. Namely, we decide that a simulation (at a given 
time) is not consistent with the observations if (i) at least 
one of the six reference quantities (cta, ctb, ARV, /a, the 
distance between the centre of mass of pop. A and pop. B 
along the line-of-sight Ay, and the distance between the 
centre of mass of pop. A and pop. B in the plane of the 
sky Axz) differs from the observed one by more than 3 
a, or (ii) Pks < 0.05 (where Pks is the probability that 
the chance deviation between the observed RV distribution 
and the simulated RV distribution is expected to be larger, 
according to a Kolmogorov-Smirnov test, last column of 
Table [2]). 


^ The intrinsic velocity dispersion for each population was de¬ 
rived as the standard deviation from the line-of-sight velocity, 
i.e. — {v))^/{N — 1) (where Vi are the line-of-sight 

velocities of each star, (v) is the average line-of-sight velocity 
and N is the number of particles). To derive d, we consider only 
the centre-of-mass motion of binary systems. 


4 
















Mapelli et al.: N-body modelling of the Gamma Velorum cluster 



Fig. 2. From top to bottom and from left to right: distribution of the radial velocities (RVs) of stars in run 1, 2, 4, 
5, 7, 11, 13, 15 and 16 at t = 4.8 Myr (see Tabled]). The simulated RVs include randomly generated observational 
uncertainties. Solid red line: pop. A (cluster 1); dotted black line: pop. B (cluster 2); dashed blue line: sum of the two 
components. We define RV as the velocity along the y axis (i.e. we assume that the line-of-sight coincides with the 
direction that maximizes t he velocity shift betw een the two clusters). Green shaded histogram: observed RV distribution 
(the same as in Figure 6 of iJeffries et al.l[^014[) . 


Run 4 (where Qa = 3) and run 7 (where Qb = 12.5) do 
not satisfy the first criterion (because of the large values of 
(TA and (Tb, respectively). While Pks is < 0.05 in runs 11 
(An = 3.5 km s“^) and 16 (/bin = 0.75). From Figure [2| 
and Table dj it is apparent that some of our models match 
the kinematics of the Gamma Velorum cluster quite well. 
Furthermore, we are able to put constraints on the most 
relevant initial parameters. In particular, if the initial ve¬ 
locity shift is An > 3 km s“^ (runs 10 and 11), then ARV 
is larger than the observed one by more than 1 cr, for the 


entire simulation. The tidal forces exerted during the in¬ 
teraction between pop. A and pop. B are not sufficient to 
quench this velocity significantly. If An >3.5 km s“^, Pks 
is very low (< 0.05). 

More interestingly, if the simulated pop. B is close to 
virial equilibrium (Qb A 2, runs 5 and 6), then its velocity 
dispersion (< 0.7 km s“^) is much smaller than the ob¬ 
served one (1.60 ± 0.37 km s“^), and the RV distributions 
of the two populations appear almost completely separated 
(Figured]). On the other hand, if Qb > 12.5, then the simu- 
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lated velocity dispersion of pop. B is too large with respect 
to the observed one (ctb ^ 2.8 km s”^ for Qb = 12.5 in 
run 7, i.e. 3 a higher than the observed value). Thus, the 
model can reproduce the observed RV distribution only if 
pop. B is significantly hotter than a virial system, but not 
too hot (2 < Qb < 12.5). This implies that pop. B is un¬ 
bound, likely as a consequence of the evaporation of the 
parent molecular gas. 

Other constraints can be put on the virial ratio of 
pop. A: if Ma 500 Mq and Qa > 3.0, then the velocity 
dispersion of the simulated pop. A 1 km s“^, run 4) is 
too large to be consistent with the observed one (0.34±0.16 
km s“^). In contrast, values of Qa < 2.0 are consistent with 
the data, indicating that pop. A is approximately in virial 
equilibrium. 

The binary fraction /bin is very important for the 
comparison with observations. In fact, unresolved bina¬ 
ries tend to broaden the RV distribution (we recall that 
in Figure [5] we plot even the simulated binaries as unre¬ 
solved, for a better comparison with the data). In most 
simulat ions, we assume /bin = 0.46 for analogy with the re¬ 
sults of iJeffries et al.l (120141 1. but we run two test cases with 
/bin = 0 (run 15) and 0.75 (run 16), respectively. The ef¬ 
fect of /bin is particularly important for pop. A. If /bin = 0 
(run 15), then the scatter in the RVs of pop. A is much 
smaller: a larger value of Qa is requested to reproduce a a if 
/bin = 0 (run 15), but the simulations still fail to reproduce 
the high-velocity tail of the RV distribution. In contrast, if 
/bin = 0.75 (run 16), the RV distribution is broader, and a 
lower value of Qa is requested to match the data. 

The characteristic radii of the two populations are more 
important for the spatial distribution than for the RV. On 
the other hand, the characteristic radii have also an impact 
on the RV, and there is partial degeneracy with the effect 
of the virial ratio. For example, if pop. B is initially spread 
over a larger radius (rn = 2 pc, run 18), for the same initial 
mass Mb, then its velocity dispersion will be lower. Thus, 
the observations can be matched only for a larger value of 
Qb, but this affects the disruption timescale of pop. B. 

Finally, the mass of the two populations is another im¬ 
portant quantity, and there is partial degeneracy with the 
virial ratio. For example, run 12 has the same parameters 
as run 1, but the mass of the two clusters is lower by a 
factor of 5. The RV distributions of pop A. and pop. B 
partially overlap in run 1, while they are well separated in 
run 12. The result is that run 12 does not match the data 
(especially ctb and ARV, Tabled]). In case of a lower initial 
total mass, larger values of Qa and Qb are requested, to 
obtain a best matching with the data. For example, in our 
simulations, if Ma + Mb ^ 150 Mq instead of 800 M0, 
we need Qa > 4 and Qb > 8 (as in run 13, Figure dj) to 
reproduce the Gamma Velorum data. 

The values adopted in run 13 match the current mass 
of the Gamma Velorum cluster quite well, but present an 
issue: they can hardly account for the formation of the 7^ 
Velorum WR-0 binary. In Table (TJ we report the masses 
of the members of the most massive binary in each run. 
The most massive binaries in runs 12 and 13 are well be¬ 
low 10 M0. Even if we account for stochastic fluctuations, 
the formation of a (30, 30) M0 binary in clusters with 
initial mass ^ 100 M0 is very unlikely. Thus, we can ar¬ 
gue either that the formation of 7^ Velorum is not related 
to the rest of the Gamma Velorum cluster (which is quite 
unlikely), or that the initial mass of the clusters was a fac- 
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Fig. 3. Proper motions in the plane of the sky for run 1 
at t = 4.8 Myr. In the main panel, red crosses: pop. A, 
black open circles: pop. B. In the marginal histograms, red 
(black) histogram: pop. A (pop. B). On the y axis of the 
marginal histograms: number of stars per bin, normalized 
to the total number of stars in each cluster (iV*/7Vtot)- The 
simulated proper motions shown in this Figure do not in¬ 
clude observational uncertainties. 


tor of 5 higher, and then the partial cluster dissolution 
due to gas evaporation has led to a smaller total mass. Our 
runs with virial ratio > 0.5 account for this effect of gas 
evaporation and ongoing cluster dissolution. Furthermore, 
our simulated YSCs evolve dynamically and progressively 
spread over an area larger than the one observed by the 
GES. This partially explains the difference between the 
mass estimated by the observations (within the 0.9 deg^ 
area surveyed by GES) and the total mass of the simulated 
YSCs in run 1 (see the discussion in Sec. 13.21) . 

Figure [3| shows the proper motions in the plane of 
the sky (the xz plane, according to our convention) for 
run 1. Pop. A and pop. B are not significantly decoupled, 
if we consider the proper motions. On the other hand, 
since pop. B is more supervirial than pop. A, the distri¬ 
bution of proper motions has a larger dispersion for pop. B 
(dx,B = 0.96 km s“^, CTz_b = 1-00 km s“^, where CTx.b 
and tTz^B are the intrinsic velocity dispersions of pop. B 
along the x and the z axis, respectively) than for pop. A 
(dx, A = 0.46 km s“^, CTz.a = 0.47 km s“^, where CTx,a and 
(tz, A are the intrinsic velocity dispersions of pop. A along 
the X and the 2: axis, respectively). Furthermore, if the plane 
of the sky is slightly different from the assumed one (i.e. if 
the relative velocity vector between the two YSCs is not 
perfectly aligned to the line-of-sight), we expect to observe 
a further difference between pop. A and pop. B. Gaia data 
will be crucial to test this feature. 
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Fig. 4. Normalized radial stellar density in run 1 (left- 
hand panels) and in run 13 (right-hand panels) at different 
times. The simulated cluster is assumed to be at a distance 
of 350 pc. Solid red line: pop. A; dotted black line: pop. B. 
The density n is calculated as stars per arcmin^, within 
concentric annuli. The normalization no was chosen so that 
the total area below the histogram of pop. A is one. The 
error bars account for Poisson uncertainties. 


3.2. Spatial distribution 

In the data published bv iJeffries et al.l (l2014ll . pop. A ap¬ 
pears to be marginally more concentrated than pop. B, 
and the two centroids are offset by 4.4 ± 3.0 arcmin 
(^ 0.45 ± 0.31 pc). Our simulations can naturally explain 
the radial offset and the different centroid of the two pop¬ 
ulations of the Gamma Velorum cluster: in our models, 
pop. A and pop. B come from two different sub-clusters, 
which might be born with different concentration and are 
expected to have different centroids. 

Figure 0] shows the normalized radial density of stars 
in pop. A and pop. B in runs 1 and 13, at different 
times. Pop. B is slightly less concentrated than pop. A. 
Furthermore, the concentration of pop. B diminishes with 
time, as a result of the fact that the second (less massive) 
cluster is being tidally disrupted by the first (more mas¬ 
sive) cluster (see also Figure [J). The probability that the 
chance deviation between the spatial distribution of pop. A 
and that of pop. B is expected to be larger, according to 
a Kolmogorov-Smirnov test, is Pks < 10“® and < 0.1 for 
run 1 and run 13, respectively. Thus, the two simulated pop¬ 
ulations follow different spatial distributions. In the other 
runs, we observe the same trend, even if the relative dif¬ 
ference between pop. A and pop. B changes (depending 
on the initial conditions). The simulated spatial distribu¬ 
tions (Figure 0]) are qualitatively similar to the observed 
radial spatial distribution of the Gamma Velorum cluster. 


as shown in Figure 11 of IJeffries et al.l ([20T3), even if we 
recall that these two figures cannot be directly comparec 0 |. 

Golumn 11 of Table 0] shows that there is an offset of 
0 — 1.5 pc between the centres-of-mass of the two simulated 
clusters in the xz plane (i.e. in the assumed plane of the 
sky) at 0 < t/Myr < 8 , depending on the initial conditions. 
For most runs, the offset is consistent with the observed 
one (~ 0.45 ± 0.31 pc). 

Finally, Table 0] also shows that there is an offset of 
0 — 18 pc between the centres-of-mass of the two simulated 
clusters along the y axis (i.e. along the line-of-sight) at 0 < 
t/Myr < 8. The colour-magnitude diagram of the Gamma 
Velorum cluster is consistent with an offset of 4±5 pc (along 
the line-of-sight), if we assume that pop. A and pop. B are 
coeval (or with an even larger offset, if we assume that 
pop. A is slightly older, as indicated by the Li depletion, 
IJeffries et al.l 120141 . Thus, the simulated offset along the 
line-of-sight is fairly consistent with the one indicated by 
the observations. 

Figure 0] shows the position of simulated stars belong¬ 
ing to pop. A (red crosses) and pop. B (black circles), pro¬ 
jected in the plane of the sky (the xz plane of the simu¬ 
lation), at three different times {t = 2.5, 4.8 and 6.9 Myr 
since the beginning of the simulation) in run 1. It is ap¬ 
parent that pop. B expands faster than pop. A, because of 
its supervirial condition. We note that stars of pop. B can 
be found at > 10 pc from the centre of pop. A. Recently, 
ISacco et al.l (|2015ll found that 15 stars, in the direction of 
the NGC 2547 cluster, i.e. about two degrees ('^ 10 pc) 
south of 7 ^ Velorum, have the same properties as the pop. B 
members of the Gamma Veloru m cluster. To com pare our 
simulations with the results of ISacco et al.l (|2015ll . we se¬ 
lect a box in our run 1 (at t = 4.8 Myr) with approximately 
the same size and the same location (with re spect to the 
centre of the cluster) as the area analyzed bv ISacco et al.l 
(1201511 . We count 16 pop. B objects (6 of which are bina¬ 
ries) in this box. Th i s is in excellent agreement with the 
result of ISacco et al.l (|2015ll . who find 15 pop. B members 
in the held of NGC 2547. This result supports our scenario, 
and strengthens the evidence that pop. B is strongly super¬ 
virial. Figure [5] also shows that our simulated star clusters 
are more extended than the area reported by Jeffries et al. 
(2014, ^ 0.9 deg^). It is reasonable to expect that the outer 
rim of the Gamm a Velorum cluster is much larger than the 
area observed bv IJeffries et al.l (120141 1. 

3.3. Collision or line-of-sight superposition? 

Our models are able to reproduce most properties of the 
Gamma Velorum cluster, provided that pop. B is signih- 
cantly supervirial. However, we may wonder whether we re¬ 
ally need a physical collision between the two sub-clusters. 
Is it possible to explain the properties of the Gamma 
Velorum cluster with a mere line-of-sight superposition be¬ 
tween pop. A and pop. B, without requiring any interac¬ 
tions between the two sub-clusters? 

® Figure 11 of IJeffries et all (120141 1 and our Figure 0] cannot 
be directly compared, since we know the intrinsic membership 
of each single star to pop. A ( i.e. cluster 1) o r pop . B (i.e. clus¬ 
ter 2), while the selection bv IJeffries et all (120141 1 is based on 
a probabilistic approach. As we have discussed in Section 13.11 
a fraction of genuine members of our simulated pop. A would 
be attributed to po p. B, based on the statistical approach by 
IJeffries et all (120141 1. 


7 

























































































Mapelli et al.: N-body modelling of the Gamma Velorum cluster 


Table 2. Main results of the iV—body simulations. 


Run t tTA CTB ARV /a o-a ARV Ay Axz Pks 


(Myr) (kins (km s (km s (km s (km s (km s (pc) (pc) 


Run 1 

2.7 

0.46 

1.93 

1.95 

0.35 

0.54 

1.29 

2.48 

1.77 

0.01 

0.72 

Run 1 

4.8 

0.43 

1.84 

1.93 

0.37 

0.50 

1.25 

2.42 

7.22 

0.16 

0.84 

Run 1 

6.9 

0.38 

1.85 

1.82 

0.35 

0.49 

1.23 

2.41 

12.37 

0.55 

0.84 

Run 2 

2.7 

0.49 

1.21 

2.66 

0.55 

0.57 

1.37 

2.51 

1.41 

0.45 

0.86 

Run 2 

4.8 

0.42 

1.66 

2.00 

0.40 

0.53 

1.33 

2.45 

6.90 

0.12 

0.84 

Run 2 

6.9 

0.47 

1.68 

2.02 

0.42 

0.50 

1.32 

2.42 

12.00 

0.06 

0.84 

Run 3 

2.6 

0.57 

1.98 

1.63 

0.32 

0.66 

1.23 

2.46 

1.92 

0.12 

0.53 

Run 3 

4.8 

0.50 

1.90 

1.67 

0.31 

0.61 

1.21 

2.41 

6.97 

0.58 

0.56 

Run 3 

6.9 

0.51 

1.85 

1.67 

0.33 

0.59 

1.20 

2.39 

12.09 

0.82 

0.68 

Run 4 

2.7 

1.04 (*) 

1.38 

2.66 

0.53 

0.95 

1.30 

2.55 

1.71 

0.06 

0.09 

Run 4 

4.8 

0.95 (*) 

1.33 

2.58 

0.51 

0.92 

1.27 

2.51 

6.94 

0.06 

0.11 

Run 4 

6.9 

0.94 (*) 

1.30 

2.57 

0.51 

0.91 

1.26 

2.50 

12.51 

0.46 

0.12 

Run 5 

2.7 

0.56 

0.55 

2.70 

0.58 

0.52 

0.73 

2.58 

1.91 

0.04 

0.24 

Run 5 

4.8 

0.51 

0.59 

2.32 

0.55 

0.51 

0.69 

2.24 

7.08 

0.25 

0.54 

Run 5 

6.9 

0.50 

0.51 

2.22 

0.56 

0.49 

0.68 

2.15 

11.54 

1.36 

0.49 

Run 6 

2.7 

0.58 

0.69 

2.74 

0.56 

0.54 

0.73 

2.65 

1.80 

0.05 

0.27 

Run 6 

4.8 

0.52 

0.67 

2.57 

0.56 

0.50 

0.69 

2.50 

7.35 

0.10 

0.37 

Run 6 

6.9 

0.51 

0.61 

2.51 

0.56 

0.55 

0.69 

2.46 

12.60 

0.16 

0.40 

Run 7 

2.7 

0.53 

2.79 (*) 

2.21 

0.42 

0.55 

2.20 

2.80 

1.62 

0.03 

0.06 

Run 7 

4.8 

0.47 

2.78 (*) 

2.21 

0.43 

0.50 

2.18 

2.80 

7.01 

0.10 

0.09 

Run 7 

6.9 

0.44 

2.79 (*) 

2.11 

0.42 

0.53 

2.18 

2.78 

11.43 

0.72 

0.13 

Run 8 

2.7 

0.49 

1.65 

1.25 

0.36 

0.55 

1.27 

1.67 

-0.50 

0.01 

0.72 

Run 8 

4.8 

0.47 

1.62 

1.20 

0.38 

0.50 

1.21 

1.64 

3.27 

0.19 

0.70 

Run 8 

6.9 

0.43 

1.57 

1.18 

0.38 

0.46 

1.19 

1.61 

6.39 

0.5 

0.40 

Run 9 

2.7 

0.46 

1.80 

1.62 

0.34 

0.55 

1.28 

2.12 

0.77 

0.01 

0.79 

Run 9 

4.8 

0.44 

1.70 

1.62 

0.36 

0.50 

1.22 

2.07 

5.48 

0.16 

0.76 

Run 9 

6.9 

0.41 

1.72 

1.53 

0.37 

0.46 

1.21 

2.04 

9.57 

0.53 

0.84 

Run 10 

2.7 

0.58 

1.35 

3.04 

0.52 

0.54 

1.33 

3.01 

3.23 

0.03 

0.23 

Run 10 

4.8 

0.55 

1.27 

2.99 

0.54 

0.50 

1.30 

2.95 

9.77 

0.15 

0.35 

Run 10 

6.9 

0.49 

1.26 

2.93 

0.52 

0.48 

1.28 

2.93 

15.58 

0.50 

0.33 

Run 11 

2.7 

0.57 

1.37 

3.49 

0.52 

0.55 

1.39 

3.44 

4.46 

0.02 

0.02 

Run 11 

4.8 

0.53 

1.37 

3.41 

0.52 

0.50 

1.36 

3.39 

11.96 

0.13 

0.04 

Run 11 

6.9 

0.48 

1.35 

3.35 

0.52 

0.52 

1.34 

3.37 

18.44 

0.70 

0.04 

Run 12 

2.4 

0.31 

0.64 

2.40 

0.54 

0.29 

0.69 

2.27 

0.87 

0.21 

0.72 

Run 12 

4.8 

0.23 

0.64 

2.40 

0.53 

0.24 

0.66 

2.25 

6.48 

0.34 

0.77 

Run 12 

7.2 

0.26 

0.72 

2.36 

0.54 

0.24 

0.70 

2.22 

11.96 

0.53 

0.81 

Run 13 

2.5 

0.33 

1.69 

1.86 

0.37 

0.48 

0.88 

2.30 

0.64 

0.01 

0.88 

Run 13 

4.9 

0.37 

1.61 

1.95 

0.39 

0.46 

0.87 

2.27 

6.31 

0.02 

0.92 

Run 13 

7.4 

0.48 

0.98 

2.35 

0.56 

0.45 

0.86 

2.26 

11.94 

0.03 

0.89 

Run 14 

2.7 

0.48 

1.95 

1.79 

0.35 

0.56 

1.35 

2.42 

-3.13 

0.02 

0.75 

Run 14 

4.8 

0.46 

1.83 

2.04 

0.40 

0.50 

1.30 

2.42 

2.30 

0.19 

0.78 

Run 14 

6.9 

0.42 

1.80 

1.96 

0.39 

0.48 

1.28 

2.40 

7.02 

0.45 

0.81 

Run 15 

2.6 

0.53 

1.07 

2.83 

0.59 

0.55 

1.38 

2.44 

1.86 

0.01 

0.42 

Run 15 

4.7 

0.40 

1.23 

2.43 

0.56 

0.46 

1.34 

2.38 

7.37 

0.01 

0.25 

Run 15 

6.8 

0.35 

1.28 

2.34 

0.54 

0.43 

1.33 

2.36 

12.78 

0.02 

0.24 

Run 16 

2.7 

0.61 

1.29 

2.87 

0.50 

0.50 

1.20 

2.75 

1.94 

0.07 

0.03 

Run 16 

4.8 

0.39 

2.19 

2.01 

0.28 

0.44 

1.17 

2.68 

6.11 

0.41 

0.03 

Run 16 

6.9 

0.37 

2.13 

1.86 

0.28 

0.42 

1.16 

2.61 

10.50 

0.44 

0.05 

Run 17 

2.7 

0.64 

1.20 

2.61 

0.57 

0.56 

1.35 

2.42 

2.12 

0.07 

0.64 

Run 17 

4.8 

0.37 

1.84 

1.83 

0.35 

0.50 

1.30 

2.42 

7.46 

0.24 

0.67 

Run 17 

7.0 

0.35 

1.75 

1.96 

0.38 

0.48 

1.28 

2.40 

12.63 

0.29 

0.67 

Run 18 

2.7 

0.57 

0.92 

2.51 

0.54 

0.54 

0.97 

2.36 

1.71 

0.02 

0.51 

Run IS 

4.8 

0.50 

0.92 

2.36 

0.53 

0.49 

0.90 

2.31 

7.11 

0.17 

0.57 

Run 18 

6.9 

0.48 

0.93 

2.33 

0.53 

0.46 

0.90 

2.26 

11.84 

0.51 

0.78 

Run 19 

2.4 

0.51 

2.11 

1.62 

0.39 

0.59 

1.58 

1.72 

1.77 

0.01 

0.76 

Run 19 

4.8 

0.44 

1.95 

1.71 

0.42 

0.49 

1.53 

1.69 

8.50 

0.17 

0.86 

Run 19 

6.8 

0.39 

1.84 

1.79 

0.44 

0.44 

1.51 

1.67 

13.70 

0.80 

0.82 

Run 20 

2.5 

0.48 

1.83 

1.81 

0.30 

0.56 

1.34 

2.20 

1.44 

0.03 

0.45 

Run 20 

5.0 

0.44 

1.72 

1.83 

0.33 

0.49 

1.29 

2.17 

8.03 

0.12 

0.60 

Run 20 

7.1 

0.39 

1.82 

1.67 

0.29 

0.49 

1.29 

2.14 

12.44 

1.00 

0.71 

Run 21 

2.5 

0.49 

1.84 

1.53 

0.34 

0.56 

1.33 

2.13 

10.40 

0.02 

0.84 

Run 21 

5.0 

0.49 

1.82 

1.62 

0.38 

0.52 

1.33 

2.07 

15.10 

0.19 

0.81 

Run 21 

7.1 

0.41 

1.84 

1.46 

0.35 

0.49 

1.33 

2.06 

19.37 (») 

0.49 

0.84 

Data^ 

- 

0.34 ± 0.16 

1.60 ± 0.37 

2.15 ± 0.48 

0.48 ± 0.11 

- 

- 

- 

4 ±5^ 

0.45 ± 0.31 

- 


Run (column 1): identifying name of the run; t (column 2): time elapsed since the beginning of the simulation; a a (column 3): 
velocity dispersion of pop. A as derived from fitting with two Gaussian components; ctb( column 4): velocity dispersion of pop. B 
as derived from fitting with two Gaussian components; ARV (column 5): RV difference between pop. A and pop. B as derived 
from fitting with two Gaussian components; /a (column 6): fraction of stars belonging to pop. A, as derived from htting with two 
Gaussian components; a a (column 7): intrinsic velocity dispersion of pop. A; as (column 8): intrinsic velocity dispersion of 
pop. B; ARV (column 9): difference between the average RV of pop. A and pop. B; Ay (column 10): distance between the centre 
of mass of pop. A and the centre-of-mass of pop. B along the line-of-sight {y axis). A negative value of Ay means that the first 
collision between cluster 1 and cluster 2 has not happened yet; Axz (column 11): distance between the centre of mass of pop. A 
and the centre-of-mass of pop. B in the plane of the sky {xz plane); Pks (column 12): probability that the chance deviation 
between the observed RV distribution and the simulated RV distribution is expected to be larger, according to a 
Kolmogorov-Smirnov test. Pks is very small (< 0.05) only for run 11 (An = 3.5 km s“^) and run 16 (/bin = 0.75). 

The values reported in the last line come from the analysis of iJeffries et al.l (120141 1. '’The observed pop. A in Gamma Velorum 
is closer by Ay = 4 ± 5 pc than pop. B, if the observed difference in the colour-magnitude diagram il.Teffries et al .1120141 is due to 
a distance difference rather than to an age difference. 

Quantities differing by > 3 cr with respect to the observed value are marked by (*). 











Mapelli et al.: N-body modelling of the Gamma Velorum cluster 


Run 1 £=2.7 Myr £=4.8 Myr £=6.9 Myr 


12 

6 

3 0 

N 

-6 

-12 


Fig. 5. Position of simulated stars belonging to pop. A (red crosses) and pop. B (black open circles), projected in the 
plane of the sky (the xz plane of the simulation), at three different times {t = 2.5, 4.8 and 6.9 Myr since the beginning 
of the simulation) in run 1. The centre of the frame is the centre of mass of the entire cluster (pop. A+ pop. B). We note 
that pop. B expands much faster than pop. A given its supervirial condition. Stars belonging to pop. B can be found at 
> 10 pc from the centre of pop. A. 



-12 -6 0 6 12 -12 -6 0 6 12 -12 -6 0 6 12 

X (pc) X (pc) X (pc) 


To answer this question, we have performed run 21, in 
which the two sub-clusters form close to each other (5 pc) 
but increase their separation rather than collide. The initial 
properties of run 21 are the same as those of run 1 (which 
matches the properties of the Gamma Velorum cluster par¬ 
ticularly well), apart from the sign of the relative velocity. 
The main kinematic properties of run 21 (i.e. the values of 
a A, ctb, and ARV) are very similar to those of run 1, indi¬ 
cating that the velocity field of the two populations is not 
severely affected by the collision. The main reason is that 
the relative velocity of the two sub-clusters is large with re¬ 
spect to their internal kinematics. On the other hand, the 
line-of-sight distance between the two sub-clusters becomes 
too large (> 19 pc, 3cr larger than the observed value) at 
t ^ 7.1 Myr. 

The same conclusion can be reached on the basis of 
a back-of-the-envelope calculation. The GES data indicate 
that there is an offset ARV = 2.15 ± 0.48 km s“^ between 
the RV of pop. A and pop. B. Furthermore, pop. A might 
be closer to us by 4 ± 5 pc. If we assume that pop. A and 
pop. B formed in the same place and then drifted away from 
each other, pop. B is now at a distance dAB from pop. A: 


dAB 


10 pc 


/ ARV \ 
V2kms-i ) 


tB \ 

5 Myr J ’ 


( 1 ) 


where Ib is the age of pop. B. Thus, even if pop. A and 
pop. B formed in the same place 5 Myr ago, dAB is already 
1 a larger than the observed displacement between pop. A 
and pop. B. If pop. B formed farther from us than pop. A, 
then we would expect a larger observed displacement be¬ 
tween pop. A and pop. B. On the other hand, if pop. A 
and pop. B formed in the same place, they should have in¬ 
teracted between each other before drifting away. For the 
same reason, if pop. B formed closer to us, it should have 
interacted with pop. A before drifting away. This favours 
the scenario of a gravitational encounter between the two 


sub-clusters. This argument is weakened by the fact that 
our measurement of the line-of-sight displacement between 
pop. A and pop. B relies on the colour-magnitude diagram 
and is quite inaccurate. Parallax measurements by Gaia will 
give invaluable hints on this point. 

Do the simulations show any other quantitative differ¬ 
ence between a simple line-of-sight superposition and a 
gravitational interaction between pop. A and pop. B? A 
physical collision leaves an imprint on the structure of the 
two sub-clusters, and especially of the lighter one. In our 
simulations, pop. B is tidally perturbed and stretched by 
the gravitational encounter with the more massive pop. A. 
The tidal perturbation can be quantified by comparing the 
half-mass radius projected in the plane where the tidal force 
is maximum (rh(j/z) or rii{xy)), with the half-mass radius 
projected in the plane where the tidal force is minimum 
{r\i{xz)). Figure ini shows that rh(yz) is significantly larger 
than rii^(xz), indicating that pop. B is severely stretched 
by the interaction. Again, this might be checked with forth¬ 
coming parallax measurements by Gaia. 

In summary, a collision scenario and a mere line-of-sight 
superposition show no significant differences in the velocity 
dispersions and in the RV distribution. The only way to 
distinguish between these two scenarios is to measure the 
line-of-sight displacement between the two populations and 
any possible tidal deformation of pop. B. 

3.4. Relaxing the assumption of line-of-sight superposition 

It is very unlikely that the two sub-structures of a YSG are 
observed exactly along the line-of-sight. How much can we 
relax the assumption of line-of-sight superimposition and 
still match the observed data? To check this hypothesis, we 
focus on our fiducial run (run 1), and change the line-of- 
sight by rotating the yz plane about the cc—axis by an angle 
0. We stop rotating when at least one of our reference quan- 
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Fig. 6. Half-mass radius of pop. B as a function of time in 
run 1. rii{xz): half-mass radius projected in the xz plane. 
Thiyz): half-mass radius projected in the yz plane (where 
the tidal force exerted by pop. A is maximum). 


tities for the comparison with the data (cta, 0-3, ARV, /a, 
Ay and Axz) differs from the best observed value by more 
than 3 cr. The hrst parameter for which the discrepancy be¬ 
comes > 3(7 is the difference between the centroids (Axz). 
At an angle 0 = 12° between the original y axis and the 
new y axis (y'), we obtain a a = 0.48 km s“^, ctb = 1.64 
km s“^, ARV=1.99 km s“^, /a = 0.40, Ay' = 7.07 pc 
and Axz' = 1.50 pc (the best value for the observations is 
Axz = 0.45 ± 0.31 pc, see Table 2). Thus, we can conclude 
that configurations in which the two sub-clusters are super¬ 
imposed or slightly offset (by an angle 9 < 12°) reasonably 
match the data. Larger angles are possible, but only for dif¬ 
ferent collision times. Figure [7] shows the RV distribution 
(along y'), the distribution of proper motions along z', and 
the projected positions (in the xy' and xz' planes) if run 1 
is rotated hy 6 = 12° about the x axis. 

The probability of observing two sub-clusters offset by 
an angle 6 < 12° is ~ 0.022. On the other hand, we recall 
that while the centroid of pop. A is quite well constrained, 
the centroid of pop. B is more uncertain, especially if pop. B 
is much more extended (as suggested bv ISacco et al.lf201^ . 
The uncertainty on the centroid of pop. B might further 
relax the request of partial line-of-sight superposition. 

4. Conclusions 

Several hints indicate that YSCs form from the merger 
of sub-clusters, born from the same molecular cloud. In 
this paper, we have simulated the collision between two 
sub-clusters. We have shown that the collision product re¬ 
produces some interesting features of the Gamma Velorum 
cluster, recen t ly ob served by the GES. The GES data 
((Jeffries et al.l 1201411 indicate that the Gamma Velorum 
cluster is composed of two kinematically decoupled pop¬ 
ulations: pop. A and pop. B. The two populations have a 
marginally different radial concentration and slightly offset 
centroids (by ~ 0.5 pc). 
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Fig. 7. Properties of run 1 (at t = 4.8 Myr) if rotated by 
an angle 9 = 12° about the x axis. Top left-hand panel: 
RV distribution (symbols are the same as in Fig. [5]); top 
right-hand panel: distribution of proper motions along z' 
(symbols are the same as in the marginal histograms of 
Fig. [21 in particular pop. A and pop. B correspond to the 
red and black histogram, respectively); bottom left-hand 
panel: positions in the xy' plane (symbols are the same as 
in Fig.lH in particular pop. A and pop. B correspond to red 
crosses and black open circles, respectively); bottom right- 
hand panel: positions in the xz' plane (symbols are the same 
as in Fig. O in particular pop. A and pop. B correspond to 
red crosses and black open circles, respectively). 


Our simulations can naturally explain the RV offset be¬ 
tween the two populations of the Gamma Velorum cluster, 
as well as their different intrinsic velocity dispersions (see 
Figure O. Our simulations also account for the different 
concentration and the different centroid of the two popu¬ 
lations (Figures [4] and [5]). The GES data suggest that 
pop. A is older by 1 — 2 Myr, based on the Li depletion. 
Our simulations are consistent with a small age difference 
between the two populations. In the scenario of a collision 
between two sub-clusters, the two populations can either be 
perfectly coeval or have formed in two slightly different star 
formation episodes. The latter hypothesis is supported by 
observations of similar age differences in other star forming 
regions (even if the question of age spread in st ar forming 
regions is still debated, see [Jeffries et al.l l201ll and refer¬ 
ences therein). 

We predict that the dispersion of the distribution of 
proper motions is broader for pop. B than for pop. A, even 
if the relative velocity vector between the two sub-clusters 
is aligned with the line-of-sight, as a result of the fact that 
pop. B is supervirial (Figure [3]). Furthermore, if a compo¬ 
nent of the relative velocity vector is normal to the line-of- 
sight, we expect a difference in the average proper motions 
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of pop. A and pop. B (e.g. Figure El). This will soon be 
tested with Gaia data. 

Furthermore, our simulations can help us reducing the 
allowed parameter space for the initial conditions of the 
Gamma Velorum cluster. For example, models with a low 
binary fraction (/bin 0) can hardly account for the ob¬ 
served RV data. Similarly, models in which the initial rela¬ 
tive velocity between the clusters is ARV> 3 km s“^ are in 
disagreement with the observed RVs. Star cluster models 
with an initial total mass <150 Mq can reproduce most 
of the observed features, but cannot account for the forma¬ 
tion of the massive 7 ^ Velorum binary system. In contrast, 
models with a larger mass (> 500 Mq) can host massive 
binary systems, similar to the 7 ^ Velorum binary. 

Finally, the simulations suggest that pop. B is not 
in virial equilibrium: supervirial models (with virial ratio 
2 < Qb < 12.5) are in better agreement with the GES 
data. The virial ratio Q is an essential ingredient of our 
models, and strongly affects the results: it is impossible to 
match the observed RV distribution of the Gamma Velorum 
cluster, without allowing Qb to assume a value >> 0.5. 
This result might indicate that the gas from the parent 
molecular cloud evaporated very fast, leaving the stellar 
population out of virial equilibrium. Thus, pop. B is about 
to dissolve in the field. Furthermore, the physical mean¬ 
ing of Q is connected with the star formation efficiency in 
the cloud. A larger value of Q corresponds to a lower star 
formation efficiency, and consequently to a stronger effect 
of gas expulsion (e.g. iTutu kovllli^Sl: iHillsIllQSfll : iLada et H] 
11984 iBaumgardt fc KrouDall2007r) . Thus, we can infer that 
pop. A formed with higher star-formation efficiency, man¬ 
aged to survive gas evaporation, and preserved virial equi¬ 
librium, while pop. B was formed in a less efficient star for¬ 
mation episode, and is not going to survive gas evaporation. 
A possible scenario is that pop. B formed from a less dense 
molecular cloud core, where star formation efficiency was 
intrinsically lower. Another possible scenario is that pop. B 
formed later than p op. A (as indicated by the Li depletion, 
iJeffries et al.l 12014) . and that the gas surrounding pop. B 
was evaporated by the radiation field of pop. A, quenching 
the star formation episode of pop. B quite abruptly. 

Since it is supervirial, pop. B expands faster than 
pop. A. Simulated stars belonging to pop. B lie at > 10 
pc from the centre of po p. A at t .^ 4 Myr since the be¬ 
ginning of the simulation. ISacco et al.l (|2015ll have recently 
found that some stars located in the NGC 2547 cluster 
have the same properties as those of pop. B. Since the 
NGG 2547 cluster lies about two degrees ('^ 10 pc) south of 
7 ^ Velorum, this result supports our scenario and strength¬ 
ens the evidence that pop. B is strongly supervirial. 

The main intrinsic limits of our simulations are the fol¬ 
lowing: we have not included the parent molecular gas (the 
effects of gas evaporation are indirectly modelled by as¬ 
suming a virial ratio > 0.5), and we do not account for the 
Galactic tidal field directly. We will refine our simulations, 
by including the missing ingredients, in a forthcoming pa¬ 
per. Our model requires a fine-tuning: we need to assume 
that most of the RV offset between the two sub-clusters 
is aligned along the line-of-sight (± 12 °; for larger angles, 
the distance between the centroids differs by > 3 cr with re¬ 
spect to the measured value). If the direction of the relative 
velocity vector between the two sub-clusters is significantly 
different from the line-of-sight, then we would expect to see 
a larger spatial offset between the two populations. If the 


Gamma Velorum cluster represents the case of a fortuitous 
alignment between velocity offset and line-of-sight, then we 
would expect to observe a number of ‘twin’ sub-clusters in 
young star formation complexes, showing a larger spatial 
offset than the Gamma Velorum cluster. The forthcoming 
data by the Gaia mission will likely shed light on this issue, 
by providing accurate proper motion and parallax measure¬ 
ments. 

Our results show that the Gamma Velorum cluster is an 
ideal test-bed to check different scenarios for the formation 
of young star clusters in the Milky Way. The RV precision 
achieved by the GES allowed us to make an unprecedented 
comparison between the kinematics of simulated and ob¬ 
served clusters. Similar kinematic data of other young star 
clusters and associations in the GES sample will be essen¬ 
tial to shed light on the formation of star clusters in the 
local Universe. 
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